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Abst  ract 


The  purpose  of  this  paper  is  to  examine  the  effect  of  the 
distortion  of  the  primary  solvation  shell  upon  the  vibrations  of 
an  ion  within  this  shell.  In  addition,  the  effect  of  an  external 
charge  (or  image  charge)  on  the  ionic  vibrations  is  also  investigated. 
An  objective  of  this  work  is  to  see  whether  it  is  possible  to  use 
predicted  changes  in  the  far  infrared  vibrational  spectrum  to  probe 
the  alterations  in  the  solvation  structure  which  take  place  when 
an  ion  is  incorporated  into  the  met al /so lut i on  interface.  In  the 
solution  phase  the  vibrations  of  the  ion  within  the  solvation  shell 
yield  a single  band,  as  was  discussed  in  the  previous  paper.  It 
is  shown  here  that  this  single  band  can  split  into  two  or  three 
new  bands  due  to  the  presence  of  distortion  of  the  solvation  shell 
and  the  effect  of  an  external  charge.  The  case  of  the  lithium 
cation  dissolved  in  dimethylsulphoxide  with  four-coordination  is 
examined.  It  is  found  that  an  external  charge  or  an  image  charge  in 
the  metal  has  a pronounced  effect  upon  the  extent  of  the  splitting  of 

the  far  infrared  band.  An  image  charge  in  the  metal,  for  example, 

- 1 ° 
yields  a splitting  of  about  16  cm  when  the  lithium  ion  is  2 A 

from  the  metal  surface.  Distortion  of  the  solvation  shell,  for 

moderate  deviations  from  uniform  spherical  symmetry,  has  a relatively 

small  effect  upon  the  system.  Nevertheless,  these  distortions 

may  have  an  important  effect  in  the  activation  process  of  the 

electron  transfer  reaction  and  other  transport  processes.  It  appears 

that  far  infrared  vibrational  spectroscopy  of  the  interface  may 

be  able  to  probe  the  structure  of  adsorbed,  solvated  ions. 


I nt  roduc t i on 


The  purpose  of  this  paper  is  to  examine  what  happens  to  the 
far  infrared  vibrational  spectrum  of  a solvated  ion  when  the 
solvation  structure  undergoes  various  kinds  of  distortion.  The 
incorporation  of  a solvated  ion  into  the  met al / sol ut ion  interface, 
for  example,  is  expected  to  be  accompanied  by  some  deformation  of  the 
primary  ionic  solvation  shell.  In  addition,  the  presence  of  an 
external  charge  or  an  image  charge  outside  the  solvation  structure 
can  have  an  effect.  It  is  therefore  of  interest  to  investigate 
these  distorted  systems  to  see  if  the  far  infrared  spectrum  might 
provide  a tool  for  the  investigation  of  interfacial  structure  and 
forces . 

The  far  infrared  vibrational  spectrum  of  the  symmetrically 
solvated  alkali  metal  cations  in  dim-nhylsulphoxide  (DMSO)  (and  in 
other  polar  solvent)  can  be  understood  in  terms  of  a simple  model. 

The  model,  which  was  investigated  in  detail  in  the  previous  paper, ^ 
consists  of  the  cation  trapped  inside  a primary,  or  inner,  solvation 
cage.  The  vibrations  of  the  ion  are  those  of  a mass  in  a spherical 
harmonic  oscillator  well.  In  the  harmonic  oscillator  limit,  it  was 
shown  that  the  classical  electrostatic,  ion-dipole  interactions  do 
not  contribute  to  the  magnitude  of  the  force  constant  as  long  as 
the  solvation  structure  retains  simple  cubic  symmetry.  Consequently, 
the  far  infrared  vibrational  spectra  provide  a probe  of  the  structure 
of  the  repulsion  forces  which  operate  in  these  systems. 

We  show  that  the  distortion  of  the  solvation  structure  can 
produce  two  or  three  new  vibrational  bands  from  the  single  band  of  the 
solution  phase  species.  In  addition,  an  external  charge  or 
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an  image  charge  exerts  a strong  influence  on  the  system. 

The  lithium  cation  in  I1M SO  is  examined.  The  primary  solvation 
cage  is  assumed  to  he  four-coordinated.  In  the  solution  phase 
it  is  assumed  to  be  tetrahedral.  At  the  met  a 1 /sol ut ion  interface 
an  image  charge  produces  most  of  the  splitting  if  the  solvation 
structure  changes  only  by  a few  hundredths  of  an  angstrom  unit. 

Two  types  of  uniform  distortion  are  examined:  (1)  uniform  distortion 
along  a C,  symmetry  axis,  and  (2)  uniform  distortion  about  a C2 
axis.  In  the  first  case,  a single  solvent  molecule  is  allowed  to 
move  in  or  out  of  the  solvation  structure  along  the  symmetry 
axis  which  passes  through  the  molecule.  The  externa]  charge  is 
located  on  the  same  C,  axis.  This  type  of  distortion  leads  to  a 
splitting  into  two  new  bands.  The  second  type  of  uniform  distortion 
allows  two  solvent  molecules  to  move  in  or  out  of  the  solvation 
structure  to  the  same  extent  along  their  tetrahedral  radii.  The 
C7  symmetry  axis  bisects  the  tetrahedral  angle  enclosed  by  the  radii 
of  the  moving  species.  It  is  also  the  axis  upon  which  the  external 
charge  is  placed.  In  this  case  the  vibrational  spectrum  of  the  original 
solution  phase  species  splits  into  three  bands  with  or  without  the 
presence  of  an  external  charge.  Two  of  the  bands,  however,  typically 
are  very  close  together  in  frequency.  Hence,  it  is  not  likely 
that  they  can  be  resolved  in  the  solution  phase  or  interfacial 
system . 

The  presence  of  an  image  charge  in  a metal  due  to  an  ion 

O 

located  2 A from  the  metal  surface  is  shown  to  yield  a splitting 
of  the  vibrational  band  of  the  order  of  16  cm  . This  result 
applies  to  a system  in  which  the  solvation  shell  remains  unchanged 
from  its  solution  phase,  uniform  structure.  As  is  shown,  changes 


* 


-4- 

in  the  solvation  shell  structure,  i.e.,  distortions,  can  increase 
or  decrease  the  extent  of  the  splitting.  Compress ional  distortions 
tend  to  decrease  the  extent  of  the  splitting.  In  these  cases 
the  ionic  vibrational  modes  appear  little  changed  from  the  solution 
phase  vibrations.  Distortions  of  expansion  tend  to  increase  the 
extent  of  the  vibrational  band  splitting.  Thus,  enhanced  vibrational 
splitting  may  be  looked  upon  as  evidence  for  a tendency  toward 
ionic  desolvation  in  the  interfacial  state. 

In  the  following  sections  these  conclusions  are  documented 
by  means  of  a detailed  analysis  of  the  distorted  ionic  system. 

Vibrations  Inside  the  Distorted  Solvation  Shell 


It  was  pointed  out  in  the  previous  paper  that  the  potential 
energy  function  for  the  solvated  ionic  system  must  be  expressed  in 
a suitable  form.^  Thus,  several  expansions  were  used  in  order  to 
express  the  small  amplitude  vibrational  excursions  of  the  ion  about 
its  equilibrium  position  in  appropriate  coordinates.  The  equilibrium 
position  of  the  ion  in  a uniform,  symmetrically  solvated  shell  is 
at  the  centre  of  symmetry.  In  this  case,  it  is  natural  to  use 
spherical  polar  coordinates  to  analyse  the  single  (s-type)  vibra- 
tional degree  of  freedom  in  the  harmonic  limit. 

On  the  other  hand,  when  there  is  distortion,  it  is  simpler 
to  use  a cartesian  coordinate  system.  The  system  of  coordinates 
is  still  appropriate  to  the  analysis  of  the  small  amplitude  vibrational 
excursions  of  the  ion  away  from  the  equilibrium  position.  However, 
the  centre  about  which  the  harmonic  oscillator  expansion  is 
constructed  need  not  necessarily  be  the  centre  of  symmetry  (if 
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one  exists),  nor  need  it  even  he  the  equilibrium  position. 

The  simplest  procedure  to  follow  when  dealing  with  a distorted 
svstem  is  to  consider  the  distortion  effects  with  reference  to  the 
symmetric,  solution  phase  system.  Therefore,  it  is  possible 
to  write  a general  (cartesian)  form  of  the  harmonic  oscillator 
potential  energy  function  as 

V = i(k  x2  + k v2  + k „z2)  + k (xy)  + k (xz)  + k (yz) 

2 xx  yyJ  zz  xy7  xz  yz  ' 

+ Xx  + Yy  + Zz  (1) 

s s s 

The  quantities  k are  the  (tensorial)  force  constants.  , Y , and 
2^  are  the  shift  forces  which  are  related  to  the  magnitude  and  direc- 
tion of  the  displacement  of  the  new  equilibrium  position  with 
reference  to  the  coordinate  origin  of  the  expansion  system  of 
coordinates.  [Throughout  this  paper,  the  coordinate  system  for  the 
expansion  in  terms  of  the  Taylor  series  is  taken  to  be  the  centre 
of  symmetry  of  the  solvated,  uniform  system.]  Thus, 

Xo  = 'Xs/kxx’  yo  = 'Ys/kyy’  z°  = * Zs/,kzz’  (2-* 

The  equations  of  motion  for  the  ion  in  the  three  cartesian 
directions  are  of  course  easy  to  determine.  They  are 


mx  = 

X 

X 

X 

-X 

1 

■ kxyy  - 

k z 
xz 

- X 

s 
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-k  x 
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- «yyy  - 
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yz 
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s 

(3) 

mz  = 

-k  x 
zx 

• kzyy  * 

k z 
zz 

' Zs 
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and 

k = k , (4) 

xy  yx  ’ 

etc.  For  the  purpose  of  determining  the  vibrational  frequencies, 
it  is  possible  to  ignore  the  , Y.,  and  Zs  quantities  in  eqn  (3). 
Thus,  the  vibrational  problem  reduces  to  one  of  finding  three 
eigenvalues  for  the  three  degrees  of  freedom.  The  equilibrium 
position  of  the  ion  in  the  distorted  system  is  given  by  eqn  (2). 

Although  the  vibrational  equations  of  motion  (3)  are  extremely 
simple  in  form,  the  major  effort  in  solving  the  problem  is  the 
determination  of  the  functional  form  of  the  force  constants  and  the 
shift  forces.  In  the  next  section  we  examine  a simple  system  for 
which  an  analytical  solution  can  be  obtained  easily.  The  system  is 
that  of  a tetrahedrallv  solvated  ion  which  is  uniformly  distorted 
by  moving  a solvent  dipole  along  a C,  symmetry  axis.  Following  this 
analysis  we  will  examine  a system  in  which  two  solvent  molecules 
are  allowed  to  move  uniformly  in  and  out  with  respect  to  a C,  axis. 

The  derivations  of  the  force  constants  and  the  shift  forces 
are  considered  in  a series  of  appendices. 


Uniform  Distortion  Aloni 


Svmmet  rv 


The  model  system  and  the  system  of  coordinates  used  are 
illustrated  in  Fig.  1.  The  solvent  molecules  are  labelled  (1) 
to  (4)  with  molecule  (4)  lying  on  the  axis  which  passes  through 
the  coordinate  origin  and  the  points  (-1,-1,1)  and  (1,1,-1).  An 
external  charge  (or  metallic  image  charge)  is  located  at  the 
point  (c,c,-c)  where  c = R-//3,  and  R.  is  the  distance  between  the 
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charge  and  the  coordinate  origin. 

for  tetrahedral  solvation,  the  sum  of  the  electrostatic  con- 
tributions to  the  potential  energy  function  from  the  dipoles  (1)  to 
(3)  is  (see  appendix  2) 

T — I -XV  ♦ X 2 + V Z ) 

1 R- 
e 

where  u is  the  magnitude  of  the  solvent  dipole  moment,  and  where  the 
location  of  the  molecule  (1),  for  example,  is  given  by  (a, a, a) 
with  a = U //?.  A single  displaced  dipole  on  the  C,  axis  is  located 
at  (-h,-b,b)  with  h = R^/ZI  where  Rj  is  the  distance  along  Cj  the 
ion  has  been  displaced.  The  ion-dipole  interaction  contribution 
from  this  term  is 


4 — l xy  - xz  - yz). 


The  sum  of  these  two  terms  is  written  as 


vdistm 


l — [ (R  /R , ) **  - 1]  (xy  - xz  - yz) 
7 » e u 


When  Rj  = Re  , Vjist(i)  = 0,  as  required. 

The  external  charge  located  at  (c,c,-c)  yields  a contribution 
of  the  form  (sec  appendix  11 

— (xy  - xz  - yz). 


Z is  the  charge  on  the  external  species.  This  term  is  added  to 

eqn  (5)  to  give  the  total  electrostatic  contribution  to  the  potential 


Next,  we  consider  the  contribution  from  the  ion-solvent 
repulsion  interactions.  The  potential  is  written  in  the  form 


p(R)  = —(1  - 2 0 / R ) exp  ( - (R  - R )/p].  (9) 

3 o2  6 

I'he  constants  B and  ; are  found  in  the  general  form  of  the  exponential 
rcpu 1 s i on  ^ 

\ r(R)  - B expf-CR  - Re)/P]  ; (10) 

B is  of  the  order  of  2.3  X 10  ^ erg  and  p is  1.9  X 10  ^ cm. 

The  distortion  contribution  is  given  by 


Vdist(r)  = + p/R)  exp  [ - ( R - Re)/p]  - (1  ♦ P/Re)/Re 


The  complete  potential  energy  function  is  expressed  as 


V = ik  (x2 


v2  ♦ z2 


+ f,(xy  - xz  - yz) 


02) 


where  f is  given  by 


c Tp 


R-0  + p/R)exp[-(R-Re)/p]  - J-  (1  + p/R  ) 


3 pe 
R 


T — [(Re/Rd)4  - 1!  + 2 e 2 / R | 


03) 


The  equations  of  motion  (3)  yield  the  following  secular 
equation  for  the  frequencies: 


3f  x + 2f3  = 0 
c c 


04) 


wi  th 


x = mu>2  - k . 


05) 


The  solutions  to  (14)  are  well  known;' 


x = 2f  , -f  , -f  , 
c ’ c c 


from  which  we  write 


w = [w?  + 2f,/.n]1/Z,  [u? 


fc/m] * ^ Itwice) 
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I-igs.  2 and  3 illustrate  the  dependence  of  the  two  frequencies 

on  the  displacement  of  the  fourth  solvent  dipole  along  the  axis. 

It  can  he  seen  that  when  the  solvent  dipole  is  displaced  a sufficient 

distance  inward,  the  vibrations  of  the  ion  a 1 ong  the  axis  are 

of  a higher  frequency  than  the  vibrations  normal  to  the  same  axis. 

for  the  ion  in  the  absence  of  an  external  charge,  this  behaviour 

changes  at  the  point  where  Rj  = Rfi . As  R^  increases,  the  vibrations 

along  the  C , axis  decrease  in  frequency  compared  to  the  vibrations 

normal  to  the  exis.  The  same  behaviour  is  seen  when  the  external 

or  image  charge  is  present  (provided  of  course  that  the  external 

charge  is  opposite  the  ionic  charge).  When  there  is  an  attraction 

between  the  ion  and  the  external  charge,  the  change-over  point 

in  the  behaviour  occurs  for  R.  < R . 

d e 

When  the  solvation  shell  is  distorted,  the  ionic  equilibrium 
position  shifts  from  the  centre  of  symmetry  of  the  solution  phase 
system  to  a new  position.  Note,  in  these  calculations  the  distance 
between  the  external  charge  and  the  centre  for  the  Taylor  series 
expansion  (to  obtain  the  harmonic  oscillator  terms)  has  been  kept 

O 

a constant  quantity.  [The  distance  between  charges  is  about  5 A 
which  corresponds  to  the  electrostatic  interaction  between  a charge 

O 

and  its  image  in  a metal  where  the  charge  is  about  2 A from  the 
metal  surface.)  The  actual  distance  between  the  ion  at  its 
instantaneous  location  and  its  image  or  the  external  charge  is 
given  inherently  in  this  work.  The  shift  forces,  the  harmonic 
terms,  and  additional  terms  in  the  Laplace  expansion  of  1/r  when 
all  taken  into  account  give  the  actual  interaction  for  the  actual 
separation.  In  constructing  the  harmonic  oscillator  system  potential 


energy  function,  it  is  possible  to  eliminate  the  linear  terms  by 
means  of  a simple  d i agona 1 i za t i on . Thus  , there  arise  constant 
terms  in  the  complete  expression  for  the  potential  energy  function 
which  together  with  the  terms  expressed  here  account  for  the 
interaction  between  charges  for  their  actual  locations.  For  the 
harmonic  oscillator  problem,  however,  these  terms  vanish  upon 
differentiation  leading  to  the  force  terms.  Hence,  they  are 
not  displayed,  nor  are  they  needed  here. 

The  shift  forces,  cf.  eqn  (2),  for  this  system  are  given  by 


X = V 


: = f 

S S 


with  (see  appendix  4) 


f = —[2^-|l  - (R  / In  . ) 3 ] + (B/p)(l  - exp  [ - ( R - R ) / p ] - e2/R- 
s /yl  R9  0 d e . 


The  new  equilibrium  positions  are 


x 0 = f s / k , y 0 = - f s / k , z 0 = f s / k . 


The  magnitude  of  the  displacement  is  given  simply  by 


R . = ZT| f . |/k 


and  t lie  direction  of  the  displacement  (along  (',1  is  given  by  the 
sign  of  f . If  f^  > 0,  the  displacement  is  away  from  the  external 
charge.  As  Fig.  4 shows,  when  the  dipole  is  displaced  inward 
toward  the  coordinate  origin,  the  ion  tends  also  to  move  in  with  it 


Uniform  Pis  to  it  ion  About  a C Symmetry  Axis 


A second  soluble  distortion  model  involves  the  uniform 
distortion  about  (instead  of  along)  a C7  symmetry  axis.  The  result 
of  this  type  of  distortion  is  the  splitting  of  the  original  single 
solution  phase  band  into  three  separate  bands.  The  analytical 
forms  for  the  frequencies  in  this  case  are  simple  in  appearance, 
and  can  he  derived  in  a straightforward  manner  from  the  system 
discussed  in  the  last  section. 

The  geometry  of  the  system  is  illustrated  in  Fig.  5.  The 
y-axis  is  chosen  as  the  C,  axis  of  reference;  obviously,  any 
axis  will  yield  the  same  results.  An  external  charge  is  located  at 


C 0 , R . , 0 ) 


The  positions  of  the  solvent  dipoles  are  seen  to  be  the 


following: 


(0,a,/7a),  «l,a,-/Ia),  (-/7a, -a, 0),  (/7a, -a, 0) 

with  a = R//7.  These  locations  can  be  generated  from  those  of 
Fig.  1 by  means  of  a rotation  about  the  y-axis  of  45°.  Thus, 


x -*•  — (x*  - z'),  y - y ' , and  z -*•  i(x'  + z')  (22) 

/I  /7 

in  the  new  frame  (henceforth  we  drop  the  prime  on  the  new  coordinates) 
The  same  transformations  can  be  applied  directly  to  the  potential 
energy  function  for  the  system  of  Fig.  1.  In  order  to  get  an 
appropriate  representation  for  this  new  system,  however,  some  care 
is  needed  to  consider  the  correct  grouping  of  terms.  Thus,  for 
example,  with  reference  to  Fig.  1,  the  ion-dipole  interactions 
for  the  dipoles  located  at  (a, a, a)  and  (-a, a, -a)  can  be  written 


, ( ih'/R"  ) ( xv  + \z  * v:  I 


4 ( m o / R " ) ( x y 


► XZ  V z 


I ho  sum  ol  those  terms  is  obviously 


3 (*: 


Similarly,  lor  the  ili  polos  located  at  (a,  a,  a)  and  ( a,  a, a)  one  gets 


3 H£  ( X 2 i . 


[f  1^!  R2  , then  Vi  *■  V.  0,  as  required  in  the  symmetric  case. 
However,  here  we  allow  R,  J R?  , and  further  applv  the  transformat 
(22)  to  the  sum  of  (23)  and  (21).  The  result  is 


5 pe 
,d  2 * 

c 


|1  (R  /R  | ) **  1 ( x 2 
1 1 1 


in  the  ( , system  o I fig.  5 . [Note,  the  transformations  (22)  are 
vibrational  principal  axis  transformations  which  have  the  effect  ol 
diagonalizing  t he  potential  energy  1 unction.) 


The  contribution  to  the  potential  from  the  external  charge 


(of  Z 


i s w r i 1 1 e n a s 


V = ( e * / 2 R ? )( x * 
e i 


finally,  hv  means  of  arguments  similar  to  those  already  used. 


the  repulsion  contribution  is  written  as 


p/Rd)exp[-(Rd-R  )/p]  - (1  ♦ p/R  ) (x 

e ‘ 


With  p(R)  defined  by  eqn  (9),  the 


2 


r 


term  is  written  as 


V 

sym 


vs  1 1 h 


k = 2p(Re)  * 2p  f Rd) 


Thus,  the  equations  of  motion  are  found  to  be 


mx  = 


-F  x = 
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k ♦ 5—  [1  - (Rp/R.)4]  * e 2 / R ? + 2f 
p 4 e cl  l r 
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-F  v 
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where  f is  defined  by 
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(Ke/Rd)4]  + e2/R?  - 2 f r z 


Jp 


i (1  ♦ p/Rd)exp[-(Rd-Re)/p]  - J-  (1  ♦ p/R  ) 
d e 


2 - z2). 
(27) 


(28) 


(29) 


(30) 


(31) 


Because  of  the  differences  in  the  expressions  for  the  force  constants 
three  vibrational  frequencies  clearly  can  be  defined: 


In  the  absence  of  an  external  charge,  and  when  R^  = Rp , the 
frequencies  (32)  are  all  equal  one  to  another.  Note,  that  in  the 
configuration  of  Fig.  5 the  solvent  molecules  between  the  ion 
and  the  external  charge  (or  image)  are  allowed  to  move.  As  the 
two  solvent  molecules  leave  the  solvation  cage  in  the  direction  of 
the  external  charge,  the  vibrations  in  the  x and  y directions  change 
There  is  a point  at  which  these  two  frequencies  are  the  same,  but 
for  which  R^  f R ; when 

3 H£[i  - (Re/Rd)4]  = - 2fr(Rd)  (33) 

Re 

the  system  of  three  bands  col  lapses  into  two. 

It  is  a simple  matter  to  reverse  the  situation  just  considered 
That  is,  it  is  easily  possible  to  maintain  the  solvent  molecules 
which  lie  between  the  ion  and  the  external  charge  in  fixed  positions 
and  to  examine  the  effect  of  the  vibrational  frequencies  as  the 
solut ion- s ide  molecules  move  further  into  the  cage  or  out  of  it. 
Retaining  the  external  charge  on  the  positive  y-axis  at  (0,Rd,0), 
but  allowing  the  solvent  molecules  at  the  positions  (-/2a, -a, 0) 
and  (/2a, -a, 0)  to  move  radially  along  their  tetrahedral  radii,  we 
get  the  following  force  constants: 

F.  . k . 3K±[1  • (Re/Rd)'I  * e»/Rj  - 2fr  - Fj 

K 

e 

Fy  ' Fy  (34) 

F.  . k . jH£[l  - (Re/Rdn  * e*/R!  * 2fr  ■ F^. 

K 


The  roles  of  and  F ^ 


are  merely  reversed  from  the  previous  case. 
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This  second  arrangement  is  more  in  line  with  the  model  discussed  in 
the  last  sect  ion . 

It  remains  to  consider  the  displacement  of  the  ion  in  this 
configuration.  for  the  situation  shown  in  Fig.  5,  the  displacement 
force  is  found  to  be 

Vs  . -i(B/D)(exp[-(Rd-Re)/pJ  • l)  * i!£  ((R./RjJ8-  1J 

e 

- e2/R? . (35) 

There  are  no  components  and  Zs,  as  symmetry  considerations 
dictate.  For  the  configuration  in  which  the  solvent  dipoles 
which  lie  between  the  ion  and  the  external  charge  are  fixed, 
it  is  easy  to  see  that 


Y ' = 
s 


-CB/p) 

I 


exp [ - (Rd-Re)/p]  - 1 


= rr[(R„ /M3  -H 


T R3 
e 


- e 2/R? 
i 


(36) 


As  in  the  last  section,  here  as  well,  yj  = - Y^/k  indicates  a 
shift  of  the  ionic  equilibrium  position  closer  to  the  external 
charge  for  a compress ional  distortion,  and  the  opposite  for  an 
expansion . 


Discussion 


The  vibrational  structure  of  a solvated  electroactive  species 
is  important  in  particular  in  the  activation  of  an  electron  transfer 


L 


reaction.  The  displacement  of  the  electroactive  species  to  a 
position  closer  to  or  farther  away  from  another  active  centre 
will  have  an  effect  on  the  rate  of  reaction.  Clearly,  deformations 
of  the  solvation  structure  which  enable  charge  transfer  species 
to  approach  closer  than  would  he  the  case  for  rigid,  symmetrical 
solvation  will  assist  the  transfer  by  allowing  the  system  greater 
overlap  potential.  Therefore,  it  is  necessary  to  have  knowledge 
of  the  vibrational  structure  of  these  solvated  systems. 

Tlie  most  likely  experimental  system  in  which  to  observe 
the  effects  of  solvation  shell  deformations  on  the  vibrational 
spectrum  appears  to  be  the  metal /solut i on  interface.  For  this 
reason,  the  proceeding  discussion  and  the  examples  have  been  given 
with  reference  to  this  system.  Nevertheless,  it  is  equally  clear 
that  distortion  effects  may  be  important  in  homogeneous  systems, 
lor  example,  the  transfer  of  an  electron  between  electroactive 
species  in  the  solution  phase  may  depend  in  part  upon  transient, 
deformed  states  of  solvation.  Such  states,  however,  are  not  likely 
to  be  observed  due  to  t he i r very  short  lifetimes.  Thus,  it  would 
seem  that  the  primary  information  about  solvation  deformation  must 
come  from  studies  of  the  inter  facial  state. 

In  this  section  several  matters  are  discussed.  Some  of  the 
discussion  is  an  elaboration  of  the  results  presented  in  the  last 
two  sections.  Some  of  the  discussion  concerns  conclusions  which 
can  be  drawn  from  the  proceeding  results.  In  particular,  the 
nature  of  the  solvation  deformation  observed  here  suggests  further 
modeling  which  can  take  into  account  the  dynamics  of  the  solvation 
shell  as  well  as  the  dynamics  of  the  ion  within  the  shell. 
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As  the  results  of  the  analyses  clearly  show,  the  vibrational 
spectrum  of  an  ion  within  a deformed  solvation  shell  can  only 
reveal  the  forces  which  operate  between  the  ion  and  the  surrounding 
solvent  dipoles.  The  analyses  presented  here  were  based  upon  the 
assumption  of  solvent  configurations.  It  certainly  cannot  be 
stated  that  one  form  of  distortion  is  necessarily  favoured  over 
another.  The  answer  to  that  question  requires  further  the  considera- 
tion of  the  interaction  between  the  solvent  dipoles,  their  images, 
the  rest  of  the  solution,  and  their  images  in  the  metal.  Such  an 
analysis  has  not  yet  been  carried  out.  However,  on  the  basis  of  the 
behaviour  of  the  ion  in  the  structures  displayed  here,  some 
tenative  conclusions  can  be  made  even  in  the  absence  of  a definitive 
interfacial  solvation  structure. 

We  note  that  unless  the  solvated  ion  approaches  very  close  to 
the  metal  surface,  there  is  no  strong  tendency  present  to  desolvate. 
The  ion  is  shifted  inside  its  solvation  cage  in  the  direction  of  the 
attractive  image.  However,  assuming  a configuration  of  close 
packing  of  solvent  molecules  on  the  metal  surface,  there  probably 
is  little  liklihood  an  ion  will  be  found  in  contact  with  the  bare 
metal  surface. 

The  expansion  of  the  solvation  structure  brought  about  by 
the  departure  (or  the  tendency  to  depart)  of  a solvent  molecule 
on  the  solution  side  of  the  interfacial  system  tends  to  take  the 
ion  with  it.  That  is,  a loose  solvation  structure  normal  to  the 
metal  surface  favours  an  increased  distance  between  the  ion  and  the 
surface.  If  the  solvent  molecule  actually  leaves,  however,  the 
ion  generally  will  remain  in  the  vicinity  of  the  remainder  of  the 
structure.  The  remaining  forces  of  attraction  between  the  ion,  its 
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image,  the  other  solvent  dipoles  are  likely  to  hold  the  ion  in  the 
interface  as  another  solvent  molecule  moves  in  from  the  solution 
to  build  up  the  solvation  structure  once  again.  On  the  other  hand, 
a compressed  solvation  structure  favours  a more  compact  system  with 
the  ion  displaced  closer  to  the  metal  surface.  It  is  not  possible, 
of  course,  to  determine  which  type  of  structure  (loose  or  compact) 
is  likely  on  the  basis  of  the  analyses  presented  in  this  paper. 
However,  in  view  of  the  possible  solvent  dipole-image  dipole 
interactions  which  are  possible,  it  is  not  unreasonable  to  expect 
that  a compact,  distorted  structure  may  be  found  in  the  interface. 

Unfortunately,  the  analysis  of  the  two  types  of  distortion 
(along  the  axis  and  about  the  C axis)  indicate  that  for 
compact  structures  the  splitting  of  the  vibrational  band  of  the 
solution  phase  system  is  not  as  great  as  it  would  be  if  there 
were  no  distortion  or  if  the  solvation  system  were  to  expand. 

For  lithium,  for  example,  the  mobile  solvent  group  would  have  to 

O 

be  compressed  0.1  A in  the  direction  of  the  metal  surface  for  the 
splitting  to  vanish.  It  is  difficult  to  say  at  this  time  whether 
such  a displacement  is  reasonable.  As  was  indicated  in  the  intro- 
duction, if  the  solvation  structure  remains  unaltered  from  its 
solution  phase  configuration,  the  splitting  is  of  the  order  of 
16  cm  A splitting  of  this  order  is  probably  experimentally 

observable . 

In  any  case,  whether  or  not  the  distortion  is  present  and 
observable  in  the  far  infrared  spectrum,  the  distortion  of  the 
solvation  structure  (if  it  exists)  will  enter  into  any  consideration 
of  the  process  of  activation  of  an  electron  transfer  reaction. 

The  dynamics  of  ionic  vibration  inside  the  solvation  cage,  as  worked 
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out  here,  apply  to  the  electron  transfer  reaction  system,  and  to 
other  solution  phase  transport  processes. 

The  examples  discussed  in  the  last  two  sections  assumed 
tetrahedral  solvation.  It  is  not  difficult  to  derive  similar 
expressions  which  apply  to  the  octahedral  solvation  case  and  to 
other  possible  solvation  structures.  Indeed,  in  the  appendices 
to  follow,  the  various  contributions  to  the  harmonic  oscillator 
limit  potential  energy  function  are  expressed  in  sufficient 
generality  to  be  applicable  to  any  distribution  of  solvent  molecules 
in  the  primary  solvation  shell. 

An  important  consequence  of  these  analyses  is  the  following. 

The  harmonic  oscillator  limit  implies  for  these  systems  that  the 

various  distortions  can  be  viewed  as  ellipsoidal  deformations. 

Thus,  in  both  cases,  i.e.,  distortion  along  C,  or  about  C2  , an 

expansion  of  the  solvation  shell  effectively  distorts  the 

spherical  system  into  a prolate  shape,  a compression  into  an 

oblate  shape.  The  dynamics  of  the  solvation  shell  definitely 

have  not  been  taken  into  consideration  here.  In  the  spirit  of  the 

Bo rn - Oppenhe imer  separation,  the  ion  has  been  allowed  to  vibrate 

within  a stationary  solvent  cage  in  any  state  of  deformation. 

In  real  systems,  clearly,  the  solvation  structure  alters  as  the 

ion  moves  through  the  solution  medium.  The  ion  probably  responds 

adiabatical ly  to  changes  in  the  solvation  structure.  The  fact 

that  the  form  of  the  distortion  is  ellipsoidal  makes  it  possible 

to  simplify  the  consideration  of  the  combined  system.  The  solvation 

2 

structure  can  be  regarded  as  that  of  a vibrating  liquid  surface. 


! 


1 

* 
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It  is  possible  to  formulate  a theory  of  the  deformation  of  this 
surface.  It  is  also  possible  to  examine  the  effect  of  the  surface 
on  the  vibrations  of  the  ion  contained  within  the  surface.  An 
adequate  treatment  of  the  solvation  surface,  and  its  vibrations, 
requires  the  consideration  of  the  primary  and  the  secondary 
solvation  shells.  Such  a treatment  is  now  in  progress.  It  suffices 
at  this  point  to  conclude  by  noting  that  once  the  solvation 
structure  and  its  dynamics  have  been  specified,  it  is  possible  to 
predict  the  nature  of  the  ionic  vibrations  inside  the  solvation 
shell  with  the  use  of  the  analysis  presented  here. 


This  work  was  supported  by  the  IJ.  S.  Office  of  Naval  Research, 
Arlington,  Virginia,  (J.S.A. 

Appendix  1:  Vibrations  in  the  Presence  of  an  F.xternal  Charge 

The  simplest  type  of  distortion  to  consider  is  that  due  to  the 
presence  of  a single  external  charge  placed  outside  the  primary 
solvation  cage.  This  charge  can  arise,  for  example,  as  an  image 
charge  in  a metal  or  it  might  be  a counter  ion  of  some  form  held 
in  place  in  a rigid  lattice. 

We  consider  an  expansion  about  the  centre  of  symmetry  for 

the  isolated,  symmetrically  solvated  ion.  Thus,  the  Laplace 
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expansion  of  the  free  space  Green  function  is 
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(1.1) 


and  R > r . V,  (r)  in  the  spherical  harmonic  function'  with  r the 
tm 

unit  vector.  Fhe  second  order  term  (P»2)  contributes  directly  to 
the  harmonic  oscillator  terms  in  the  vibrational  potential  energy 
function.  The  first  order  terms  are  associated  with  the  actual 
displacement  of  the  ion  away  from  the  centre  of  symmetry  associated 
with  the  symmetrically  sc'vated  ion.  These  first  order  terms  will 
he  considered  in  appendix  4. 

The  second  order  electrostatic  contribution  can  be  expressed 
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The  coordinates  (X,Y,Z)  locate  the  position  of  the  external  charge 
with  reference  to  the  coordinate  origin  at  the  centre  ot  symmetry 
of  the  solvated  ion.  The  location  of  the  ion  with  respect  to  the 
same  coordinate  origin  is  given  by  (x,yfz).  Upon  differentiation 
of  these  ionic  coordinates  (x,y,z),  the  contributions  to  the  force 
field  are  easily  found. 

The  simplest  case  to  consider  places  the  external  charge 


along  the  principal  axis: 


viz.,  at  R on  z . Thus  , 


As  a result,  it  is  possible  to  write 


k . = k0  + 2ZiZ2e2/R3 


(along  the  z-axis) 


(1.4a) 


k = k0  - Z i Z2e2/R' 


(normal  to  the  z-axis) 


(1.4b) 


where  the  k are  force  constants,  and  k^  is  the  force  constant  with 
reference  to  the  symmetric  system  in  the  absence  of  the  external 
charge . 

Assume,  for  example,  that  the  additional  charge  is  an  image 
charge  in  the  metal  for  a solvated  ion  in  the  interface.  Then, 


k+  = k0  - e2 / 8D* 
k = k0  + e2/16D3 


(1.5) 


where  D is  the  distance  between  the  ion  and  the  metal  surface.  If 
D is  2 A,  the  original  lithium  hand  splits  into  two  bands  separated 
by  16  cm  1 . 


Appendix  2:  Ion-Dipole  Contributions  to  the  Harmonic  Oscillator 

Potential 


It  was  shown  in  a previous  paper1  that  for  symmetric,  cubic 
solvation  there  is  no  second  order  (£=2,  r ) contribution  to  the 
potential  energy  function  arising  from  classical  electrostatic 
terms.  This  ceases  to  be  the  case  when  the  solvation  structure 
departs  from  simple  cubic  symmetry.  In  this  section,  the  general 
formula  for  the  second  order  contribution  to  the  potential  which 
arises  from  the  interaction  between  an  arbitrarily  placed  dipole  and 
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ion  located  near  the  coordinate  origin  is 


The  Car  1 son  - Kushb rooke*  expression  to  order  0(r2)  was  given 


earl i er  as 
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Uiere  a is  the  dipolar  unit  vector,  p is  the  magnitude  of  the  dipole 
V'liith  is  located  at  R = (X,Y,Z).  The  transformation  of  this 
equation  to  cartesian  coordinates  yields 
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Apporulix  3 : Contributions  from  the  Exponential  Repnision 


Given  the  Fourier  transform  of  a function,  it  is  possible  to 
express  the  second  order  Taylor  series  expansion  term  as' 
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where  ro  is  the  expansion  point.  The  Fourier  transform  has  the 
form 
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where  j n ( x ) is  the  spherical  Bessel  function  of  the  first  kind,7* 
and  ( £ i H 2m  i m2  | H 3m  3 ) is  the  Clebsch -Gordan  coefficient.'1 

The  form  of  the  repulsion  interaction  used  in  the  last  paper 
was  exponential  of  the  form' 

Vr(r)  = B exp [- (R-a)/ p]  (3.4) 


where 


R = 2epp/R* 


(3.5) 
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and  p is  a constant  of'  the  order  of  1 . 9 x 10'9  cm.  R is  the 

equilibrium  solvation  radius  for  the  symmetrically  solvated  system. 

If  we  consider  an  expansion  of  the  repulsion  about  the  centre 

ot  symmetr>  lor  the  hulk  phase  solvated  ion,  then  r0  = 0,  and  in 
1 

general , 
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wliere  k^fx)  is  the  modified  spherical  Ressel  function  of  the  third 
kind.3  For  H=0,  A = 0,  and 

f2 (0)  = ^r2 (B/p2)exp(Re/p) | (R/pJkj (R/p)  - 3kQ  CR/p ) ] 

= £r2 (B/p2) (1  - 2p/R) exp [ - (R- R ) / p ] (3.7) 

The  only  other  contribution  to  t,  comes  from  £=2,  A = 2.  Thus, 
l2(2)  (r)  = 2rrr2  (B/p2)  (1100  | 20)  2 (2200  | 00)  | (R/p) k,  (R/p) 

5k2 (R/p) ] X (22u -M I 00)V  (R)Y?  (r)  (3.8) 

Ihe  summation  in  eqn  (3.8) can  be  written  as 


r 


- 


Y 2 ) ( x 2 


(22w-u  |00)  V 2ij  (R)  Y2  _u(r) 


.•>  / b 
1 6 7T 


{(X' 


- y2) 


(3Z: 


R ) (3z 


7 , 2 


llonco  , 


r2)  ♦ 4 [ (XY) (xy ) - (XZ)(xz)  ♦ (YZ) (yz) ] }/R2 r2 

(3.9) 


t:(2)(r)  = y(B/pR’)(l  + P/R)exp( -(R-Re)/p] {(X2-Z2)x2 
* (Y 2 - Z 2 ) y2  ♦ ( 3Z 2 - R2 ) z 2 + 2 [ (XY) xy  + (XZ)xz 
+ (YZ)yz ] } (3.10) 

When  this  term  is  combined  with  t,^()^,  we  find 

= T[(p+q)x2  + (p+r)y2  + (p  + s)z2l  + 6xy  + ex z + cyz  (3.11) 

with 

P=  (B/3p2) (1 -2p/R)exp[ - (R-Re)/p]  (3.12) 

and 


q = f(X2-Y2) 
r = f (Y2 -Z2 ) 
s = f (3Z2 -R2 ) 

(3.13) 

6 * fXY 
e = fXZ 
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f = (B/pRJ) (l*p/R)exp[- (R-R  )/p] . (3.14 
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The  first  order  Taylor  series  term  for  the  repulsion  inter- 
action is  calculated  with  the  use  of  the  integral  form  of  the 
Taylor  series  as  discussed  in  ref.  1.  Thus,  with  the  Fourier 
transform  of  cqn  (3.4)  given  by 


v ( k ) = 8uB  exp  (R  /p)  5 (4.1) 

c p(l/p2+k2)2 


the  first  order  Taylor  series  term  is  given  by 


1 1 (r)  = r • | Vr\'r  (r  ,R)  ] 


= — r-fd’k  k v (k ) exp  ( - i k • ( r - R ) ] _ (4.2) 

( 2 tt ) 3 ~ j ~ ' ~ ~ 


As  was  discussed  in  ref.  1,  we  make  use  of  separate  expansions  of 
the  arguments  in  the  exponential  in  the  second  line  of  eqn  (4.2). 
The  resulting  integration  can  be  carried  out  directly.  The  result 
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t , (r)  = - i=-f  H / p) exp (R  / p)  r Y . f r ) Y . (R)  [ 3k  ] (R/ p) 

1 ' m 

- (R/p)k2(R/p)] 

= (B/pR)exp[- (R-Re)/p] (Xx  + Yy  + Zz).  (4.3) 


Note,  the  value  of  Rp  is  the  equilibrium  solvation  radius  with 
respect  to  the  symmetric  system. 

The  electrostatic  terms  are  somewhat  more  complicated  in 
appearance.  Again,  use  is  made  of  the  Car 1 son - Rushbrooke  expansion 
to  express  the  shift  terms  in  an  appropriate  form.  The  first  order 
(2=1, r)  term  contributes  to  the  shift;  it  is  given  by 


1 (CR) 


ue- 


r (4tt  ) 


3/2 


R3  3/5 


(Y1_1(r) [/5Y1.1(a)Y22(R) 


/^Y10  ( a ) Y 2 1 (R)  + Yn(a)Y20(R)]  + Yu(r)  [Y1.1(a)Y20(R) 


- /3Y10(a)Y21  (R)  + /6Yn  (a)Y2_2(R)]  + Y1Q(r) [ -ZSYj.j (a)Y21 (R) 


♦ 2Y10(a)Y2Q(R)  - /3Yu(a)Y2.1(R)]} 


(4.4) 


Rewriting  this  expression  in  terms  of  cartesian  coordinates  gives 


V1 (CR) 


j^[dxx  + dyy  + dzz 


(4.5) 


- - 
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Vv  1 t 


dx  - 1 

a (ZX2-Y2-Z2)  + 3a  XY 
x y 

+ 3 a , X Z ] / a R 2 

dy  = 1 

3a  XY  + a ( 2Y2 -X2  - Z2) 
x y 

+ 3a, YZ]/ aR2 

(4.6) 

d:  = | 

3a  XZ  + 3a  YZ  + a (2Z2 
x y z 

-X2 -Y2 ) ] /aR2 

Finally,  an  external  charge 

Ne  located  at  (X' 

, Y ' , Z ' ) ye i Ids 

X 1 (ex) 

= No'  [ X ' x + Y'v  ♦ Z 
(R'  ) 3 

' z ] - 

(4.7) 

Urns  , 

upon  combining  terms, 

we  can  write 

V1  = X 

x + Y y + Z z 
s s s 

(4.8) 

where  X , Y , and  Z stand  for  the  shift  forces;  these  forces  are 

s s s 

given  specifically  by 


Xs  = (B/pR)exp[- (R-Kc)/p]X  - ped^R3  ♦ Nc2X'/(R')3 

Ys  = (B/pR)exp[ - (R-Re)/p]Y  - pedy/R3  ♦ Ne 2 Y ' / ( R ' 1 3 (4>9) 

Zs  = (B/pR)exp[- (R-Re)/p]Z  - pcd„/R3  + Ne2Z ' / (R ’ ) 3 . 
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Figure  Captions: 

1.  The  configuration  of  a solvated  ion  with  a coordination 
number  of  4 at  the  surface  of  a metal.  The  solvent  molecules  S-^ 
to  lie  in  a plane  coplanar  with  the  surface.  The  axis 
is  normal  to  the  surface  and  passes  through  the  ion  and  the  solvent 
molecule  . D is  the  distance  between  the  ionic  Expansion  centre 
and  the  metal  surface.  The  actual  location  of  the  ion  depends 
upon  the  shift  forces  (see  text  for  details). 

2.  The  dependence  of  the  vibrational  frequencies  in  the  system 

of  Fig.  1 as  a function  of  the  distance  the  solvent  molecule  S4  is 

displaced  along  the  axis.  These  frequencies  have  been  calculated 
for  the  system  in  the  absence  of  an  image  or  external  charge.  The 
line  - - A - - is  associated  with  the  vibrations  along  the  axis. 

The  other  (two-fold  degenerate)  vibrations  are  normal  to  the 
C,  axis.  The  frequencies  are  three-fold  degenerate  at  the  solution 
phase  equilibrium  point  for  the  distance  between  the  ion  and  the 
solvent . 

I 

3.  The  dependence  of  the  vibrational  frequencies  in  the  system 

of  Fig.  1 as  a function  of  the  distance  the  solvent  molecule 

is  displaced  along  the  axis.  These  frequencies  have  been 
calculated  for  the  system  in  the  presence  of  an  image  charge 

O 

located  at  a constant  distance  from  the  metal  of  2 A.  The  frequencies 
associated  with  the  line  --A--are  those  along  the  C,  axis. 


A 
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4.  The  displacement  of  the  ion  with  reference  to  the  centre  of 
symmetry  of  the  solution  phase  system.  The  displacement  along 

the  line  --A--  is  associated  with  the  system  in  the  presence  of  an 
image  or  external  charge. 

5.  The  con f igurat ion  of  solvent  molecules  about  the  ion  in  the 
case  where  the  C7  axis  (colinear  with  the  system  y-axis)  is  normal 
to  the  metal  surface.  Two  of  the  solvent  molecules  lie  in  a plane 
which  is  coplanar  with  the  metal  surface  and  which  lies  between 
the  surface  and  the  ion.  The  other  two  solvent  molecules,  on  the 
solution  side  of  the  ion,  lie  also  in  a plane  coplanar  with  the 
surface.  As  in  Fig.  1,  D is  the  distance  between  the  ionic 
expansion  point  and  the  surface. 


6.  The  dependence  of  the  vibrations  of  the  system  of  Fig.  5 
as  a function  of  the  displacement  of  the  solvent  molecules  in 

pairs.  The  solid  line  represents  the  vibrations  in  the  y- 

direction  along  the  Cj  axis.  The  vibrations  along  the  line  --+-- 
are  along  the  x-axis  and  --A--  are  along  the  z-axis.  This  applies 
to  the  case  for  which  the  molecules  between  the  ion  and  the  surface 

< 

move.  If  the  solvent-side  molecules  are  displaced  in  pairs,  then 
the  vibrations  along  --+--  correspond  to  vibrations  in  the  z-direction 
and  along  --A--  correspond  to  vibrations  in  the  x-direction. 

7.  The  dependence  of  the  vibrations  of  the  system  of  Fig.  5 

in  the  presence  of  an  image  charge.  The  interpretation  of  the  lines 
is  the  same  as  in  Fig.  6. 
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8.  The  displacement  of  the  ion  as  a function  of  solvent  molecule 
displacement.  Curve  --A--  indicates  the  displacement  of  the  ion 
in  the  system  in  the  absence  of  an  image  charge.  The  solid  line 
indicates  the  displacement  in  the  system  in  the  presence  of  an 
image  charge.  In  both  cases,  the  solvent-side  molecules  of  Fig.  5 
are  mobile. 
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